
Title : OFDM channel estimation and tracking for multiple transmit antenna 



Description 

Field of the invention 

5 This invention relates to communication using Orthogonal Frequency Division 

Multiplexing ('OFDM') and, more particularly, to channel estimation and tracking in 
OFDM communication. 

Backoround of the invention 

This invention relates to communication using Orthogonal Frequency Division 
10 Multiplexing ('OFDM') and. more particularly, to channel estimation and tracking in 
OFDM communication. 

Background of the invention 

OFDM communication has been chosen for most of the modern high-data 
rate communication systems (Digital Audio Broadcast - DAB, Terrestrial Digital 

15 Video Broadcast - DVB-T, and Broadband Radio Access Networks - BRAN such 
as HIPERLAN/2, IEEE802.11a/g, IEEE802.15.3a. for example, and is considered 
for future wide-band telephony standards, referred to as "4G"). However, in most 
cases the receiver needs an accurate estimate of the channel impulse response. 
Also, in the context of wireless local area networks ('WLANs'), current data rates 

20 (54Mbps on top of the physical layer) are foreseen to be insufficient for very dense 
urban deployment, such as for hot spot coverage. This is the motivation for IEEE 
to propose and specify in the scope of the IEEE802.11n (the former High 
Throughput Study Group) solutions for very high data rate WLANs (targeting at 
least 100Mbps on top of the medium access control ('MAC') layer) in the 5GHz 

25 band. Another area of investigation is that of higher frequency bands where more 
spectrum is available such as the 60GHz band. 

One way of achieving higher data rates is to provide the system with multiple 
antennas both at the transmitter and at the receiver. By doing so, it is possible to 
increase the quality of the communication link by exploiting the spatial diversity 
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dimension using for instance Space Time Block Codes ('STBC'), or to increase the 
spectral efficiency of the system by transmitting simultaneously different streams 
using Spatial Division Multiplexing. Therefore, Multiple Transmit Multiple Receive 
(MTMR) antenna systems are strong candidates for next generation WLANs and 
5 certain other OFDM communication systems. 

In the so-called Cyclic Prefix OFDM (CP-OFDM) modulation scheme, each 
OFDM symbol is preceded by a guard interval that is longer than the channel 
impulse response (CIR) and a cyclic prefix or postfix, hereinafter refen^ed to 
collectively as a cyclic affix, is inserted at the transmitter in a guard interval 

10 between consecutive OFDM symbols, the cyclic affix consisting of samples 
circularly replicated from the useful OFDM symbol time domain samples. The 
cyclic affix enables very simple calculation for the equalisation at the receiver, 
where the cyclic affix is discarded and each truncated block is processed, for 
example using Fourier Transform (usually Fast Fourier Transfonn (FFT)), to 

15 convert the frequency-selective channel output into parallel flat-faded Independent 
sub-channel outputs, each coresponding to a respective sub-carrier. For 
equalisation purposes, numerous strategies exist. Following the zero forcing 
approach, for example, each sub-channel output is, unless it is zero, divided by 
the estimated channel coefficient of the corresponding sub-canier. 

20 In the Zero Padded OFDM (ZP-OFDM) modulafion scheme, as described for 

example in the article by B. Muquet. Z. Wang. G. B. Giannakis, M. de Coun/ille. 
and P. Duhamel entitled "Cyclic Prefixing or Zero Padding for Wireless Multicanier 
Transmissions" IEEE Trans, on Communications, 2002. the cyclic affix is replaced 
by null samples. This solution relying on a larger FFT demodulator, has the merit 

25 to guarantee symbol recovery irrespective of channel null locations when the 
channel is known. 

However channel estimation and tracking remains an issue, especially in the 
presence of high mobility or high frequency and data rates. Like other digital 
communication systems. OFDM modulation encounters problems at high Doppler 
30 spreads, which occur notably when the user is moving fast, for example in a car, 
or even at pedestrian speeds when investigating the area of higher frequency 
bands where more spectrum is available such as the 60GHz band. Accordingly, 



CML01 meP-EPC spec VGnal 



2S August 2003 



-3- 



the channel impulse response needs to be constantly tracked and updated, 
especially in the presence of high Doppler spreads. 

It would be desirable for the OFDM modulation system to keep all the 
advantages of classic OFDM and additionally allow very simple and (semi-)blind 
5 channel estimation at the receiver. Semi-blind channel estimation means that 
substantially no additional redundancy is added to the system with respect to 
classic CP-OFDM, and therefore no bandwidth for data transmission would be 
lost; however, semi-blind channel estimation can be realized thanks to 
deterministic sequences known at both the transmitter and the receiver sides, as 
10 long as there is no substantial bandwidth loss for data transmission. Such a 
system would be advantageous in low-mobility scenarios and would make OFDM 
systems applicable to high-mobility scenarios as well. 

Our co-pending European Patent Application EP 02 292 730.5 describes a 
communication method in which the CP-OFDM time domain redundancy is 

15 replaced by a pseudo-randomly weighted deterministic sequence which leads to 
the so called Pseudo Random Postfix OFDM (PRP-OFDM). The advantages of 
being able to use ZP-OFDM are preserved and low complexity channel estimation 
at the receiver is made possible. Note that PRP-OFDM does not impact the 
achieved useful data rate and spectral efficiency compared to the classical CP- 

20 OFDM modulator, apart possibly from transmission of small amounts of data for 
the calculation of pseudo random parameters, since the only modification is the 
affix content, thus the low complexity channel estimation possible at the receiver 
side is also semi-blind. 

Our co-pending European Patent Application describes the application of 

25 PRP-OFDM to single transmit antenna systems and It is desirable to apply 
comparable techniques to MTMR systems, capable of using more than one 
transmit and/or receive antenna. 

Summan/ of the invention 

The present invention provides a method of communication using Orthogonal 
30 Frequency Division Multiplexing, a transmitter and a receiver as described in the 
accompanying claims. 
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Brief description of the drawings 

Figure 1 is a schematic block diagrao) of a transmitter in a communication 
system in accordance with one embodiment of the invention, given by way of 
example, 

5 Figure 2 is a schematic block diagram of a receiver in the communication 

system whose transmitter is shown in Figure 1, 

Figure 3 is a diagram of signals appearing in operation of the modulator of 
Figure 2, and 

Figure 4 is a diagram of a moving average Doppler module for the 
1 0 demodulator of Figure 1 , 

Detailed description of the prefenred embodiments 

Figure 1 and Figure 2 show an OFDM communication system in accordance 
with one embodiment of the invention comprising a transmitter including an OFDM 
modulator 1 and a receiver including an OFDM demodulator 2, the transmitter and 
15 the receiver communicating over a communication channel 3. 

The OFDM communication method of this embodiment of the present 
invention enables estimation and tracking of the Multiple Input Multiple Output 
CMIMO') channels in coherent multiple transmit antenna / multiple receive antenna 
{'MTMR') systems, without any specific limit to the number of transmit (TX) and 
20 receive (RX) antennas and without imposing any particular Space-Time Code 
(STC). Data and affix vectors are independently encoded by two STC and enable 
a semi-blind estimation of all the MIMO channels exploiting only the order-one 
statistics of the received signal. 

In the following description, lower (upper) boldface symbols will be used for 
25 column vectors (matrices) sometimes with subscripts N or P emphasizing their 
sizes (for square matrices only); tilde (*-') will denote frequency domain quantities; 
argument / will be used to index blocks of symbols; « (^ will denote Hermitian 
(transpose) operations. 

Figure 1 depicts the baseband discrete-time block equivalent model of an A/- 
30 camer PRP-OFDM MTMR transceiver with Nt transmit and Nr receive antennae. 
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The communication system is described with reference to Space-Time (ST) block 
codes CSTBCs') but it will be appreciated that the invention is also applicable to 
other code systems. In the transmitter, the initial serial bit stream of constellation 
symbols ^OW),- is converted to a set of vectors In a seriaMo- 
5 parallel converter (not shown); the /th Nx1 input digital vector s(j) is then 
modulated by an Inverse Fast Fourier Transform (IFFT) matrix in a 

transfomier4. where [F^],^ =--j=wS ,0^k<N,0<.l<N and W^^e'''^^'', 

The resulting Nx1 time domain vector s{j) is processed by a suitable ST 
encoder matrix fW in an encoder 10. as shown in Figure 1. creating outputs S(i)= 

10 9A{s[iNX'-.^i}N,+N,-l))--{sXiM + k)j<l<N,fi<k i is the block 

number / is the number of the TX antenna and n^iM-^k indexes the outputs in 
Figure 1 . It will be appreciated that, at least in the context of STBCs. M can differ 
from Nt] In particular, the STBC may lead to rectangular S(i), that is to say that M> 
Nt. For the sake of simplicity, we assume in the following description that M 

15 operates on inputs (s(iAAj, - ssCw, However, the invention can be 
straightfonvardly applied to iW with other numbers of inputs. In this embodiment of 
the invention, the s^{iM+k) are lineariy preceded in a precoder 11 by a zero- 



Uf{n)^T2pSf{nll<l<N, resulting in the output vectors q/(«). The output vectors 
q,(/2) are converted to a series signal by a parallel-to-series converter 6, a pseudo 
random postfix being inserted in the signal into guard intervals between each 
consecutive OFDM symbol to produce a series digital signal ^,(«) on the fth TX 
25 atenna. The series digital signal ^,(w) is then converted to an analogue signal 
SfiO by a digital-to-analogue converter 7 and transmitted over the channel 3. 



padded OFDM ('ZP-OFDM') preceding matrix T^p, where J2p = 




and 



20 



u,(/i)=T;^s;(4;</<i\^,. 

The affix contents v,(») are added to the data symbols 



CML01 134EP-ePC spec Vfifw/ 



28 August 2003 



-6- 



More particularly, In a preferred embodiment of the Invention, the postfix that 
is added in the guard interval comprises a pre-calculated suitable vector that is 
independent of the data and that is weighted by a first factor a^, and a second 
factor Wi{kY In one embodiment of the invention, the first factor a^^ is different from 

5 one time-domain OFDM symbol to another and is known both to the transmitter 1 
and to the receiver 2, so that any time domain (cyclo-)stationarity (leading to 
strong undesired frequency contributions at the repetition frequency) is avoided. In 
another embodiment of the invention, in which the symbols are coded in blocks, 
the first factor is different from one time-domain OFDM symbol block to another 

10 but is the same for each symbol of the same block. The second factor wif/cj 
enables one of the transmit antennas to be distinguished from another. 

With an OFDM modulator in the transmitter functioning in this way, semi-blind 
channel estimation in the receiver can be done simply and at low arithmetical 
complexity. In particular, the receiver can constantly estimate and track the 

15 channel impulse response without any loss of data bandwidth compared to CP- 
OFDM, other than the transmission of PR-calculation parameters. Moreover, the 
demodulator at the receiver can have advantageous characteristics, ranging from 
very low arithmetical cost (at medium perfomnance) to les low arithmetical cost 
(with very good system perfomnance). 

20 As described in our copending European Patent Application refen^ed to 

above for the single antenna case, several choices for the first factor are 
possible. It is possible to choose of any complex value. However, any a^^with 
1 9^1 leads to performance degradation compared to prefen-ed embodiments of 
the invention. 

25 It Is possible to limit the choice of , somewhat less generally to a complex 

value with |a^| = l . This choice usually leads to good system performance, but the 
decoding process risks to be unnecessarily complex. Preferred values of the first 
and second factors are described in more detail below. 

Preferably, the first factor is pseudo-random. In one embodiment of the 

30 present invention the first factor is deterministic and is calculated both by the 
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modulator 1 and the demodulator 2 using the same algorithm and parameters that 
are stored in memory both in the transmitter and in the receiver. In another 
embodiment of the present invention, initialisation parameters for the algorithm are 
transmitted between the transmitter and the receiver to avoid systematically using 
5 the same sequence for the first factor , In yet another embodiment of the 
present invention, the first factor is communicated from the transmitter 1 to the 
receiver 2, which still represents an acceptable overhead in certain circumstances. 

In the embodiment of the present invention shown in Figure 1 and Figure 2, 
the affix is detemiinistic and is a Dxl postfix vector c treated in an encoder 12 by 
10 a specific ST encoder matrix 'W^ which outputs the Dxl vectors Cf{n\l <1<N,. The 
way in which ensures Identification of the complete MIMO channel is described 
in more detail below. The postfix vectors c,(n)L7 < / < A^, are then lineariy precoded 



adders 15: q,(w)=u,(w)+V;(»), 1<1<N, to produce signals 16 for transmission. 

The signals received at the receive antennas are the transmitted signals 
multiplied by the Channel Impulse Response ('CIR') H^^^ and with the addition of 
noise and interference n^{rt). Let be a PxP circulant matrix whose first row is 

20 given by [hj0l0,-'0,h^{L^l)r-AM^ where h,„=[hM'",hjL-l\0,'-0Y is 
the Px1 channel impulse response between the /th transmit and the /nth receive 
antennae; D is chosen such that D > L-7 . Define H^' as the lower triangular part 
of H,^ including the main diagonal which represents the Intra-Symbol-lnterference 
(ISI); Hj^' shall contain the upper triangular part of H,^ representing the Inter- 

25 Block-Interference (IBI). such that H^=H^+ H™ . Therefore, the received signal 
vector on the /nth antenna, l<m<N^ is. given by: 



in a precoder 13 by a ZP-OFDM matrix Tp, where 




, to produce 



zero-padded postfix vectors as shown at 14 in Figure 3. 

The resulting vectors V/(/j) are finally added to the data symbols u,(/i) by 
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r„W = i:(H;rqX")+H-q, 



where n^{n) is an zero-mean additive white independent identically distributed 
gaussian noise temi. 

As shown in Figure 2, The demodulator 2 at the receiver comprises an analogue- 
5 to-digital converter 7 that converts the signals r„(0 received at the receive 
antennas to digital signals, a serial-to-parallel converter 8, which converts the 
received digital signals to received vectors r„(/?), and a demodulator and equaliser 
9 that uses decoding matrices corresponding to the encoding matrices 'W^and M to 
estimate the Channel Impulse Response CIR and demodulate the OFDM signals. 

10 In the following description of the operation of the receiver, an order-one 

channel estimation algorithm is described first, assuming the channel to be static. 
Then, the effect of Doppler is introduced for the mobility case and the 
corresponding channel estimator in the Minimum Mean Square Error (MMSE) 
sense described. 

15 First the received vector r„(n) is expressed in an exploitable form for 

channel estimation. For that purpose, let Hj^ be the DxD drculant matrix of first 
row [h,Mfi-'OML'l\-Mit We define nfj^^ and such that 

H,^ =H;f'^ +H™*^. The signal r^(/i). received during the nth OFDM symbol on 
the nAh antenna, l<m<N,,\s equal to: 



where s,^(/j), s^j{n), n„„(«), n„;(») are respectively the first D and last D 
samples of s,(«) and n„(w). 

Equation 1 indicates that a superimposition of the various postfixes 
convolved by the corresponding channels is interfering with the useful data. An 
25 easy independent retrieval of each of the channels based on the sole observation 



20 r„(n) = i; 




Equation 1 
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of the postfix contributions is obtained tlirough isolation of each postfix convolved 
by its related channel. As detailed below, a way to achieve that condition is to 
perfomi a Fast Fourier Transfonn on the postfixes in the demodulator and 
equaliser 9 using a weighting ST block coding scheme W of the postfix c 
according to the following postfix generation process : 



Equation 2 

where ® is the Kronecker product and c , a{n) are respectively the detemiinlstic 
10 postfix and the pseudo-random weighting factors introduced in our co-pending 
European Patent Application EP 02 292 730.5 for the single antenna case. The 
pseudo-random weighting factors a{n) are used to convert the detemiinistic 
postfix c into a pseudo-random one. Note that a new set of detemninistic weighting 
factors is introduced, and gathered in the MxN, matrix W conresponding to the 
15 matrix W used for encoding the postfixes in the transmitter encoder, with 
[wL^.y =>v,(jt), 0<k<MJ<l<N,. W is used to remove the interference 
between all transmitted postfixes and thus is invertible in this embodiment of the 
present invention: full column rank (rank(w)=: N, ). In the following description, we 
choose W orthogonal for this embodiment of the present invention, such that 
20 W''W = I^. 

With the assumption of a static channel, an order-one channel estimator in 
the demodulator and equaliser 9 functions as follows. The first and last D samples 
of r„{n) are denoted respectively by r^^oi^) and r„;(«). By setting n-iM-^k and 
assuming the transmitted time domain signal s,(«) to be zero mean for all /. we 
25 use Equations 1 and 2 to compute for each /c, 0 ^ * < Af , the following Dx1 vector: 
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r„,{iM + k)+r„,,{iM + k + l) 



Equation 3 



Next, d^=E[d*(i)] is defined as the expectation of d*(z). Due to the 
deterministic nature of the postfixes, it can be verified from Equation 1 that: 



< =ilH,r +H2'^k(*>=iH>,(*)c 



Equation 4 



^1 

z 

5 Thus the MDxl vector d« =[(d^f , -,(^^"0^1^ can be expressed for each 

receive antenna as: 



3, 



H^w,(0)c 








= (W®lJ 




H,>,(il/-7>: 







Equation 5 



Since W is chosen orthogonal, multiplying each d„, l<m<N, by 
(W®!^)'' in the demodulator and equaliser 9 removes completely the 
10 interference between channel contributions Hj^ , ]<^1<N,. 

Once the interference between channel contributions is removed the 
estimation algorithms of the single-antenna case of our co-pending European 
Patent Application EP 02 292 730.5 can be applied in the demodulator and 
equaliser 9: 

15 Hj2.c = = FifC^F^h^ Equation 6 

where is a DxD circulant matrix with the first row [c(o),c(D-/).--,c(/)], 
=diag{F^c}. and hf„ represents the D first coefficients of h,„. Hence, the 
estimate hj^ of the time domain channel impulse response hf, in the demodulator 
and equaliser 9 is obtained by multiplying H^, by F^Q'F^. 2<l<N,J^m^N,, 
20 Note that C'J is a diagonal matrix that is known to both the transmitter and 
receiver and can thus be precalculated. Subsequently, hf^ is preferably 

transformed to the Px1 frequency domain vector = F4l^,o;;^fhf..This MIMO 
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channel estimation (i.e. estimation of all channels between any transmit and any 
receive antenna) is used to space-time decode and equalise the received data 
signals, as described in more detail in examples below, such that the transmitted 
data signals can be recovered. 

5 The above channel estimator can be extended to further improve reception in 

mobile environments by using any Doppler model, and by minimizing any 
perfonnance criterion. An example is now given, in a preferred embodiment of the 
present invention based on the introduction of a Doppler model; the estimator alms 
at minimizing the Mean Square Error (MSE). 

10 The Doppler module shown in Figure 4 is introduced in the demodulator and 

equaliser 9 to modify the order-one autoregressive model for the Channel Impulse 
Response ('CIR') between transmit antenna / and receive antenna m separately: 
hf,(«) = J<,(275roAr)h^(/2-7)+h^(7i) where y^O is the 0th order Bessel function, 
is the Doppler frequency. AT is the MTMR PRP-OFDM block duration and 

15 h^(w) is zero-mean complex Gaussian of constant variance. The same CIR 
correlations as the ones provided by the known Jakes model are obtained even in 
the presence of large Doppler frequencies. This is achieved by forcing the 

correlation E[(h^ h,^ {k)\= J, [inf^ {k - n)AT)\/k >n. This way the 
approximation Jo{WD{k-^)^T)^{Jf^{2Trfg,AT)f^''\ inherent to the order-one 
20 autoregressive estimation, is avoided. This modification leads to the following 
moving average estimation: 

h/^(«) = Z^o(2»c/z>AAr)hj2:(/i-t) Equation 7 

As for the order-one autoregressive model, so-called process-noise vectors 
h^(w) are introduced assuming '^i{n%i{n^\-^fi for ni^k and hjS,(o)=hj2.(o) 
25 being the CIR to be estimated. Assuming that CIR hf;(») is estimated based on Z 
noisy observations (/) = {d,^ (it), / < ifc < i + z}. the expression 

d,^(i)=[(W®I^)''d„(/)J^,(/~7)D<it</Z) results from the convolution of the /th 
block postfix by channel h^(w) conrupted both by themnal noise and the OFDM 
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10 



15 



data symbols. The OFDM data symtxils are assumed zero-mean and independent 

of same variance as the postfix samples. d„(i)=\[dl,(i)J , - Mm''(j)jJ contains 
the received symbols after equalization of the pseudo random weighting factor of 
the postfixes (Equation 3). Thus, d,„({) can be expressed as follows: 



5 d^(0=(J®C^] 



+ n, where 



Equation 8 



J = 



0 

1 



... 0 
'•. 0 



J„{2irf^{Z-l)AT) J,{2iir„iZ-2)AT) - 1 



where n gathers the thermal noise and the interference from the OFDM data 
symbols. In order to guarantee the unit variance of each channel realisation, the 
norm of hi^(«)t/i=0, -.Z-i is chosen such that: 



SWt EfteWf h£,(n)l = 1. A: = 0,....Z-1. Equation 9 

#f=0 L 



It can thus be verified that the optimum estimator of h^(o)=h^(o) In the MMSE 
sense is given from Equation 8 by: 

h/m(0 = M[(j®C^)DR^(j®C^)'' +R„l''d„(i) Equation 10 

with: 

M = [ir •,^(2;5^^(Z-l)Ar)]®(Rs^(o)C^ 

and R. is the auto-correlation matrix of the vector a . 

Since in practice the channel power profile is usually not known, in that case 
the assumption is made that Rjjj^^^j for all n. The real gain g„ is 

introduced for respecting the power constraints of Equation 9. 
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The above description presents generic channel estimation in tlie 
demodulator and equaliser 9 in accordance with embodiments of the present 
invention for both relatively static and high mobility environments. Their use for two 
STBC systems will now be described. 



includes modulators using pseudo-random postfixes at the transmitter and also 
equalizer structures derived for the MTMR case from those described for the 
Single Transmit Single Receive (STSR) case in our co-pending European Patent 
Application EP 02 292 730.5 based on the transformation of the received PRP- 

10 OFDM vector to the ZP-OFDM case. The system is described for the case of 
Nj=2 transmit and = 1 receive antennas, although it will be appreciated that 
the system is applicable to other numbers of antennas. The ST encoder operates 
over NfXM vectors with = A/ = 2 . Since A^^ = 1 , the subscript l<m<N/\s not 
used in the following analysis. Perfect knowledge of the channels h^, l^/<iV, Is 

15 assumed but it vAW be appreciated that the system is capable of woricing with 
imperfect channel knowledge. 

At the transmitter, a 2x1 ZP-ST encoder fw is used, which takes two 
consecutive OFDM symbols s(2/) and s(2i+l) tofomi the following coded matrix: 



20 where the permutation matrices Pjf are such that, for a Jxl 



Since the channel has been estimated, as for the single antenna case 
described in our co-pending European Patent Application, it is always possible to 
retrieve the MTMR ZP-OFDM signals from Equation 1 by subtracting from the 
25 received signal the known PRP contribution: 



5 



The first embodiment of STBC is based on ZP-OFDM decoding. The system 




Equation 11 



vectora = [a(0), -,a(J-l)7, wehave ^JJa}^ =a((7-/7+n)mody). 



r^W = rW-X[H;^'v;(/7-l)+Hf v,W] Equation 12 
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which leads to r^(n)= 2=,H/TzpS,(w). Note that i) no constraint has to be set on 

Wfor the symbol recovery, ii) the PRP interference cancellation procedure 
proposed is generic and can be applied to any suitable ST encoder M. 

A suitable detection algorithm Is applied to the signal described by Equation 
5 12 by the demodulator and equaliser 9. Noticing that P^Tzp =Ta,PjJ, we denote 
by D, =<liag^,}, D, =diag^,}, n(20 = Fpn(2i) and n(2/ + l)=F^P^n*(2/+l); then 
if we sv\^ch to the frequency domain by computing y(2i)=F;,r^(2i) and 

y(2/ +l)=Fp fi^^r^ili+ljf . exploiting the fact that H, = F^D,¥p . 1 < / < 2 , we can 
write: 



10 



y(2/) l^fD, 5, T FpT^s(2/) 1[ n(2/) I 
y(2/+l)J [d; -D;JLFpTzps(2/ + 1)J [n{2i+l)\ 



Equation 13 



Y 

= D 



where D is an orthogonal channel matrix. Thus multiplying ^inf y(2i+lfrby 
achieves the separation of the transmitted signals s(2i) and s(2/+l), and it 
can be shown that fiill transmit diversity is achieved. Note that the separation of 
15 signals allows the same equalisation schemes to be used in this embodiment of 
the present invention as in the singie-antenna case described in our co-pending 
European Patent Application EP 02 292 730.5. 

The second embodiment of STBC system is based on equalization of the full 
received block by diagonalisation of pseudo-circulant channel matrices. The ST 

20 data encoder £W used in the demodulator and equaliser 9 is based on a version of 
the single antenna system described in our co-pending European Patent 
Application EP 02 292 730.5 modified to enable the equalization structure that is 
detailed below and outputs blocks of N^xM vectors with =M = 2. 5W and 
are specified such that they generate the following matrix 

25 Q(0 = {q; (2/ + 1 < / ^ 2, 0 < * < 2} at the antenna outputs: 
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Q(0= 




■s*(2/ + l) T 
a(2/ + l)c*J 
' s{2i) 1 
a(2/ + l)c'J 



Equation 14 



fp being a permutation matrix defined as previously (inversing the order of the 
vector elements), a{i) is complex with |a(/]| = I being pseudo-random complex 
weighting factors as defined in our co-pending European Patent Application for the 
5 single antenna case with a{2i+\) = fi^{i)a{2i), and j3{i) = a{2i-2)/a{2i). Q^(^) is 
defined as: 



The Dxl postfix vector c is chosen such that it has Hemnitian symmetry, 
that is to say that the complex conjugate of the vector read backwards is equal to 
10 the original c . As in our co-pending European Patent Application, the channels are 
represented by PxP pseudo-circulant channel matrices Hf^'^ I^/^2. These 
are Identical to standard circulant convolution matrices with the upper triangular 
part multiplied by the scalar factor /3{i) . in other wonds uf'^ = uf + . 




With R(/)=[r^(2/) r^(2/+l)f and the noise matrix 
1 5 N(/) = [n^ (2i) (2/ + l)f , the received signals over Af = 2 symbol times are given 



as follows: 



X[Hrq,(2,-l)+H;%(2/)] 
R(0= +N(/) 
X&irq,{2i)+H;"q,(2/+1)] 



With R(i) = [f^(2/) f''(2i+l)f , the following operations are perfomjed at the 



demodulator and equaliser 9 on the received vectors: 
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cf)(i) = --[2SRH2i)}H(^^+23W20}Hf Jlpc. 

in this embodiment of STBC system, tlie data symbols are separated in the 
5 demodulator and equaliser 9 by premultiplication of R(/) by the Hemnitian of the 
upper channel matrix W^(j): 




with H^^'^=Hf^'^(Hf''>)^+Hf'^(Hf'^)^. The equalisation based on pseudo 
circulant channel matrices is then perfonned as presented in our co-pending 
10 European Patent Application for the single channel case. The PRP-OFDM postfix 
based blind channel estimation is performed based on r(/) as presented above. 
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